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Abstract 

Interface free energy is the contribution to the free energy of a system due 
to the presence of an interface separating two coexisting phases at equilibrium. It 
is also called interfacial free energy or surface tension. The content of the paper is 

1) the definition of the interface free energy from first principles of statistical me- 
chanics; 

2) a detailed exposition of its basic properties. 

We consider lattice models with short range interactions, like the Ising model. A 
nice feature of lattice models is that the interface free energy is anisotropic so that 
some results are pertinent to the case of a crystal in equilibrium with its vapor. 
The results of section [2] hold in full generality. 

1 Interface free energy in Statistical mechanical 
1.1 Definition of the interface free energy 

Consider a physical system at equilibrium in a vessel V, at a first order phase transition 
point, where two bulk phases, say A and B, coexist (for example, an Ising model at zero 
magnetic field and low temperature). If, when we bring into contact the phases A and 
B, the state of the system is inhomogeneous and there is spatial separation of the two 
phases, then at the common boundary of the two phases emerges a spatially localized 
structure, called the interface. In spite of the low dimensionality of these interfaces and 
their negligible contribution towards the global overall properties of the physical system 
their presence is essential for a wealth of important processes in physics, chemistry and 
biology. Here we consider only systems at equilibrium, which is a rather severe restriction. 

The interface free energy is a thermodynamical quantity and it is best explained 
when we consider the macroscopic scale, in which the length of the vessel containing 
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the system is the reference length. Under this scale the interface is well-defined and 
localized. In the case of a flat interface perpendicular to a unit vector n it is described 
mathematically by a plane perpendicular to n, which separates the bulk phases, and 
the state of the system is specified above this plane by giving the value of the order- 
parameter of one of the bulk phases, say A, and below the interface that of the other 
bulk phase. The interface free energy or surface tension (per unit area) r(n) describes the 
thermodynamical properties of the interface at equilibrium. How does one obtain r(n) 
once the interatomic interactions of the system are given? We can answer this question 
so that we get interesting information about T(n) only for few models. However, the way 
of defining r(n) is quite general and can be applied in principle to most systems, and its 
origin can be traced back to the monumental work of J.W. Gibbs, On the Equilibrium 
of Heterogeneous Substances (1875-1878). In statistical mechanics the thermodynamical 
functions are obtained by computing the partition function. The system is enclosed in 
a vessel V; taking into account the interations of the system with the walls of V we 
can write an expression for the overall free energy of the system. The basic postulate is 
that we can separate the various contributions to the overall free energy F(V) at inverse 
temperature (3 into two parts, up to a small correction term; one part is proportional 
to the volume of V, which is the bulk free energy of the system, and another one is 
proportional to the area of the surface of V, which is interpreted as the wall free energy. 
Thus, at a point of first order phase transition, when only phase A is present, 



where Z A (V) denotes the partition function of the system for phase A, \V\ the volume 
of V and \dV\ the area of the boundary dV of the vessel. A similar expression holds 
for phase B. The bulk terms /buik(^4) and /buik(-S) are the same because the system 
is at a first order phase transition point, but the surface terms /„ a ii(^4) and f wa i\(B) 
may be different. Under specific conditions on the walls, we can obtain (macroscopic) 
inhomogeneous states with planar interfaces separating the two coexisting bulk phases. 
In such cases there is an additional contribution to the overall free energy and we postulate 
that the free energy can be written as 



Fab(V) = ~\nZ AB (V) = f hulk (AB)\V\+ f wan (AB)\dV\ + T(n)\I(n)\ + o(\dV\) , (1.2) 



The term |/(n)| = 0(|<9U|) is the area of the interface perpendicular to the unit vector 
n, and a is the proportion of the walls in contact with phase A. Since the system is 
at a first order transition point fbuik(AB) = fbuik{A) = fbuik(B). Extracting wall free 
energies is not easy, but it is not necessary to do this if our postulate is correct, because 
we can eliminate the terms involving f W3l \{AB) and f wa ii(AB) by considering the ratio of 
partition functions, 



f a{V) = ~lnZ A (V) = f hulk (A)\V\ + U a n(A)\dV\ + o(\dV\) , 



(1.1) 



with 



f wall {AB) = af mll (A) + (1 - a)f wall (B) . 



1 Z AB (V) 
p n Z A (V)«Z B (V) 



—a 



r(n)|/(n)| + (|<9V|). 



(1.3) 



Notice that HI .31) is always a term of order OddVl). 
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An obvious difficulty in getting r(n) is that we must know the values of thermody- 
namical parameters of the system for which there is phase coexistence. Indeed, for other 
values of these parameters the system has only one bulk phase and there is no interface. 
Hence the surface tension is non-zero only for a specific range of values of the thermody- 
namical parameters of the system. This is why in many situations one proceeds differently 
in Physics. One models directly the interface in order to bypass these problems and then 
the interface free energy is simply identified with the free energy of the model for which 
one has standard methods for evaluating it. This is often an adequate way to proceed, 
but it cannot be applied always. For example when one is studying how the coexisting 
phases are spatially distributed inside the vessel V, we cannot avoid considering the free 
energy of interfaces between coexisting phases. 



1.2 A paradigm, the Ising model 

We implement the ideas of section 11.11 for the Ising model, for which the mathematical 
results are the most complete. We consider the three-dimensional Ising model. The 
two-dimensional case is also of interest. Let Z 3 := {t = (ti,t 2 ,t 3 ) : tj e Z} and 

A LM := {t E Z 3 : max(|ti|, |t 2 |) < L, \t 3 \ < M} . 

For each t G A^m, we introduce a variable o~(t), called hereafter spin, which specifies the 
state of the system at t; there are only two possible states labeled by 1 and —1, hence 
o~(t) = ±1. The configurations of the system are written a and are specified by giving 
for each site t the value of o~(t). The spin o~(t) interacts with an external real magnetic 
field h, the interaction energy being —ha(t), and it interacts with the spin o~(t'), the 
interaction energy being —J(t,t') a{t)a{t'). The energy of the system is equal to 

H LM (a) = -l E J(t,t')a(t)a(t')- £ ha(t) . 

te\ L M t'eA LM t&\ LM 

We choose the coupling constants J{t,t') = 0, except if t,t' are nearest neighbors, in 
which case J(t, t') — J > 0, a fixed value. This interaction favors the alignment of spins, 
since the energy is minimal when a(t) = o~(t'). To model the interaction of the system 
with the walls, we introduce an inhomogeneous magnetic field J'rj(t) acting only on the 
spins located at the boundary of the box A^m, 

W2 M (a) :=- J'vitHt) , 

teA LM -. \t 3 \=M 

or max(|ti|, | i 2 1 ) = L 

with J' > and T](t) = 1 or T](t) = — 1, but fixed. Different kind of walls are modeled by 
specifying different values for T](t) (see below) and choosing different values for J'. The 
overall energy of the system is Hlm + Wlm- 

According to statistical mechanics, the state of the system at equilibrium and at 
inverse temperature f3 is the Gibbs measure, so that 

Prob(a) = ^ " where Z\ M = re^^feO) . 

LM „, 
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The normalization constant Z r [ M is the partition function and the overall free energy of 
the system in the box A LM is 

1 



F2 M ((3,h,J') :=--\nZl M . 



At the thermodynamical limit (infinite- volume limit), the bulk free energy per spin 
is independent on the choice of J' > and of 77, 

W/?, h) = Km {2L \ 1)d F 2 L (^ h, J') . 

The model exhibits a first order phase transition with the magnetization as order-parameter 
if the external magnetic field h = and the inverse temperature f3 > /3 c (d), where f3 c (d) 
is the inverse critical temperature of the d- dimensional Ising model (0 < /3 C (3) < (3 C (2)). 
For h = and (5 > f3 c {d) the spin-flip symmetry of H LM is broken. There is a phase 
with positive magnetization m*(j3) and another with negative magnetization — m*((3)] 
the bulk free energy /buik(/3, h) is not differentiable at h — 0, 

< m*(j3) = ^/buik(/3,/i)U=o+ = -^/buik(/3,/i)U=o- • 

If we choose i](t) = +1, respectively 77 (£) = —1, that is pure boundary conditions, then 
one can prove that for any value of J' > the interactions with the walls favor the bulk 
phase with positive, respectively negative, magnetization. We are also interested in the 
case of mixed boundary conditions, which is related to the emergence of a planar interface 
perpendicular to n. Let n = (ni,n 2 ,n 3 ). We set 



-1 if tini + t 2 n 2 + t 3 n 3 > 0, 
■1 if tini + t 2 n 2 + t 3 n 3 < 0. 



Thus i] n (t) = 1 iff t is above or in the plane 7r(n) perpendicular to n and passing through 
the origin, otherwise r] n (t) = —1. We set Z£ M '■— Z V LM . 



+++++++++++++++++ 




The field J'i] n acting on the boundary spins of A LjL . 
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1.2.1 Macroscopic limit and definition of r in the Ising model 

We suppose that the external magnetic field h = and (3 > (3 c (d), d = 2,3. We choose 
< a < 1 and we partition the set An into cubic cells Cj of linear size L a ; for each cell 
Ci we introduce the empirical magnetization (averaged magnetization over the cell Cj) 

mcM ■= |Cr X 5>(t). 

teOi 

We scale all lengths by L~ l , so that the distance between neighboring spins is L" 1 . 
In the fixed reference macroscopic box V = {(xi,X2,x 3 ) G R d : \xi\ < 1} we make a 
coarse-grained description of the state of the system by specifying the coarse-grained 
magnetization profile, 

Pl(x\o) := m c (a) if (Lx x , Lx 2 , Lx 3 ) G C i: V (xi,x 2 ,x 3 ) G V. 

By definition the probability of the profile Pl{x\<I_) is the joint probability of the block- 
spins mcXgi) induced by the Gibbs measure. The macroscopic limit is the limit when 
L^ 1 tends to 0. If we choose rj = +1, when L~ l tends to 0, then for any value of 
J' > the probability measure on the density profiles becomes concentrated on the unique 
magnetization profile p(x) = m*(/?); this constant profile describes the macroscopic state 
of the +- phase of the model. The macroscopic limit corresponds to the regime of the law 
of large numbers in probability theory. 

To obtain r(n) we use (11.31) and the fact that by symmetry of the model Z^ L = Z^ L , 
so that we only need to compare Z^ L and Z" L . We set 

1 1 Z n 

r(n) :=-——- — — lim — —r In • (1.4) 
v ; /3|vr(n) HV\ i-ootf-i Z+ L v ; 

The term |7r(n) fl V\ is the area of the intersection of the plane vr(n) with V. One can 
prove: the limit (11.41) is independent on J' > J, and for [3 > [3 c {d) the function r(n) 
verifies properties a), b) and c) of {\2.7\i ; in the macroscopic limit the measure on the 
density profiles is concentrated on the unique magnetization profile 



Pn{X) 



+m*((3) if x is above vr(n) 
-m*(j3) if x is below vr(n). 



This profile describes a macroscopic state with a planar interface perpendicular to n. 
Therefore r(n) can be interpreted as the free energy of that interface perpendicular to n. 

The condition J' > J is important, because for some values of J' < J and j3 the 
physics near the walls of the system is different: a surface phase transition may take 
place and portions of the interface may be pinned to the walls. As a consequence of this 
phenomenon in the macroscopic limit the interaction of the system with the walls given by 
r] n may not induce an interface perpendicular to n. For example, in the two-dimensional 
case, the macroscopic state may have an interface making an angle with the vertical walls 
of the vessel, whose value is given by the Young-Herring equation, so that (II. 4p may not 
be equal to r(n), or, if J' is small enough and the macroscopic box is a square, then the 
whole interface may even be pinned to the walls, so that there is no interface through the 
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macroscopic box. In such cases the limit (11.41) depends on J'. The condition J' > J has a 
simple physical interpretation; it ensures that the walls of the box V are in the complete 
wetting regime, so that the interface cannot be pinned to the walls. In the literature 
the standard choice for ferromagnetic models is J' = J, so that (11.41) gives the correct 
definition of r(n). These results illustrate the fact that it is important to choose correctly 
the interactions of the system with the walls in order to use definition (11.31) . One must 
avoid the possibility of pinning of the interface to the walls. Any wall interactions which 
induce a macroscopic state with an interface perpendicular to n and such that otherwise 
( 11. 31) is independent of the chosen interactions are admissible for defining the interface 
free energy. 

Several other definitions for r(n) have been proposed for the Ising or similar models. 
Most of them involve a ratio of partition functions and are based on the same pattern 
leading to 1 1 1.4ft so that we shall not review them here (see section [3] for references). A 
possibility of avoiding the above problem with the walls is to suppress (partially) the 
walls of the system by taking (partial) periodic boundary conditions. Then one imposes 
a condition implying the existence of a single planar interface perpendicular to n. There 
are also variants of ( 11.41) where one considers a box Alm instead of All and take first 
the limit M — > oo before taking L —>■ oo. When J' < J this limit may give a different 
answer as the limit (11.41) . On the other hand, if J' > J, then one can take the limits 
in any order, first L —>■ oo and then M — > oo or vice-versa, or simultaneously L —>■ oo 
and M — > oo. The reason is that the walls are in the complete wetting regime and the 
interface is not pinned to the walls. In general it is not easy to show that reasonable 
definitions give the same value for r(n). 

The surface tension for the two-dimensional Ising model can be computed exactly. 
Onsager computed the interface free energy for n = (0, 1), 

/?r((0, 1)) = 2(K -K*), /3> P c (2) and r((0, 1)) = , otherwise, 

where K* is defined by exp(— 2K*) = t&nhK and K = (3 J. Onsager did not use the 
definition ( 11.41) ; the computation of r((0, 1)) defined by ( 11. 4ft is due to Abraham and 
Martin-L6f. The full interface free energy has been computed by McCoy and Wu. 

1.2.2 Inequalities for r in the Ising model 

We set in this subsection 

n(0) := (0, - sin 0, cos 0) and r(0) := r(n(0)) . 
There are two inequalities which relate the interface free energy and the order-parameter, 

r(0)<2(m*(/3)) 2 , (1.5) 

and 

M>W. (1-6) 

These inequalities indicate that r(0) > iff there is a phase transition. Since r(0) is a 
minimal value, we have r(n) > iff h = and (3 > (3 c {d). 
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We now introduce the step free energy r stcp . This quantity is interesting only for 
d = 3. Let T]*(t) be defined as 



By definition, if d = 3, 



-1 if * 3 > or if * 3 = and * 2 > 0, 
-1 if t 3 < or if t 3 = and t 2 < 0. 



1 - 1 - zi L 



^ := -^i^2LTI ln zt with n* = (0,0,1). 

One has r stcp > and we state two important inequalities. The first inequality 

r(0)-<r(O) > |sinfl|r step (1.7) 

gives information about the non-differentiability of r at = 0. Positivity of r s t ep implies 
that t{6) is not differentiable at 9 = since r has a minimum at 9 = and 

The physical consequence of the non-differentiability of r is explained at the end of 
section 12.2.31 non-differentiability of r implies the existence of a facet for the equilibrium 
shape W T defined in (12.11) . The second inequality relates the step free energy of the 
three-dimensional model to the two-dimensional interface free energy, which we write 
here r 2 (9), 

rstcp > r 2 (0) . (1.8) 

Therefore, if (3 > (3 C {2) then T step > 0. 



2 Basic properties of the interface free energy 

We examine in this section what are the basic properties of r and we discuss the ther- 
modynamical stability of interfaces. 



2.1 Convexity of the interface free energy 

We assume that r(n) > for each unit vector n. We consider the three-dimensional case. 
Elements of the Euclidean space E 3 are written x = {x\, x 2 , £3); we denote the Euclidean 
scalar product by 

( x|y ) := x x y\ + x 2 y 2 + x 3 y 3 , 

and the Euclidean norm by ||x||. 

By convention r(n), with ||n|| = 1, is the physical value of the interface free energy 
of an interface perpendicular to n. It is convenient to extend the definition of r to E3, 
positively homogeneous function, by setting 

r(x) := ||x||r(x/||x||). 
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We introduce the half-space 

H(n) := {x: (x|n) < r(n)}, 

whose boundary is the plane dH(n) = {x: (x|n) = r(n)}. Notice that H(n) = H(tn) 
for all t > 0. Besides we define the equilibrium shape W T as the convex set, which is 
the intersection of the half-spaces H(n), that is 

W T := {x: (x|n) < r(n), Vn}. (2.1) 

The next argument is due to Herring |H (1951)] ; it is reproduced almost in its original 
form. Let T(Ao, Ai, A2, A3) be the tetrahedron with vertices Aq, A%, A2, A3. The face 
opposite to the vertex Aj is denoted by A^ it contains all vertices Aj, j ^ i, and its area 
is |Aj|. Let rij be the outward unit normal to the face Aj. We have 

|A |n + |Ai|ni + |A 2 |n 2 + |A 3 |n 3 = . 

Set 

|Ai| |A 2 | |A 3 | 
n := -n = Tr-j-ni + 7-7-7*12 + 77-^3 ■ 
|A)| |^o| |^o| 

We compare the free energies |A |T(n) and |A 1 |r(n 1 ) + |A 2 |r(n 2 ) + |A 3 |r(n 3 ), or which 
is the same, r(n) and 

( |A 2 | |A 3 | 
r r(ni) + — — -r(n 2 ) + — — -r(n 3 ) . 



lAnl An IA 



0| 



Let rri!, m 2 , m 3 be reciprocal vectors to rii, n 2 , n 3 , defined by the relations ( m^n.,- ) = 
if i 7^ j and ( mj|rij ) = 1 otherwise. We have ( rrij|n) = |Aj|/|A |, so that 

3 |A-| 3 
J2jX~\ T ^ = (5Z r ( n m il n ) = (y| n ) • 

i=i ' °l i=i 

The vector y = Y^=i T ( n i) m i verifies the identities { y | ) = r(rij), % = 1,2,3, i.e. y is 
the intersection point of the three planes dH(xii), <9/f(n 2 ) and <9if(n 3 ). We have three 
cases. 

(1) . If ( y|n ) < r(n), then the intersection of the plane dH(n) with W T is empty, since W T 
is a subset of if (ni)nif (n 2 )flif (n 3 ). A (hypothetical) planar interface perpendicular to n 
with interface free energy r(n) can be deformed using the tetrahedron T(Aq, Ax, A2, A3) 
by lowering its free energy. Hence, at nonzero temperature such a planar interface with 
interface free energy r(n) is unstable and cannot exist at equilibrium since the free energy 
is minimal. Therefore we must have 

|A |r(n) < |Ai|r(ni) + |A 2 |r(n 2 ) + |A 3 |r(n 3 ) . (2.2) 

(2) . The case (y|n) = r(n) is a borderline case. Here the intersection of dH(n) with 
W T is at most one point. If dH(n) n W T = 0, then one also show that the interface 
is also unstable and cannot exist at equilibrium. If dH(n) PI W T = {y}, then y 6 
W T fl OH (ni) n dH(j\2) n dH(n 3 ) is called a corner of W T . In this case dH(ja) is a 
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support plane of W T at y, that is W T C H(n) and W T fl dH(n) 3 y, but it is not a 
tangent plane of W T . Notice that if W T has no corner, then at equilibrium we only have 
the next case (3). 

(3). If (y|n) > r(n), then 

|A |r(n) < \&x\T{n x ) + |A 2 |r(n 2 ) + |A 3 |r(n 3 ) . (2.3) 

We can repeat the same argument with a right prism whose base is a triangle with 
vertices A ,Ai,A 2 and whose length is very large. Let be the side of the triangle 
opposite to the vertex A4, \£i\ be its length, and rij the outward unit normal to the side 
£i in the plane of the triangle. Again we have 

|4|n + |£i|n 1 + |£ 2 |n 2 = 0, 

and we set n := — n = |£i|/|4>l n i + |^2|/|^o| n 2- In the plane spanned by ni and n 2 let 
rri! and m 2 be reciprocal vectors to rix and n 2 . Then 

2 \£-\ 2 
Yl TF~\ T ^ = ( r(n i )m i |n ) = ( z|n) . 

The vector z = ^ 2 =1 r(rij)mj belongs to the intersection of the planes dH{ja.x) and 
dH(n 2 ). We can draw similar conclusions as above. If (z|n) < r(n), then the interface 
dH(n) with interface free energy r(n) is unstable and does not exist at equilibrium. 
Hence, at equilibrium, 

\£o\r(n) < Ki|r( ni ) + \£ 2 \r(n 2 ) . (2.4) 

If ( z|n) = r(n), then either dH(n) n W T = or dH(n) n W r ^ 0; in the latter case W T 
has an edge or a corner and dH(n) is a support plane (but not a tangent plane of W T ). 
If (z|n) > r(n), then 

|4|r(n) < |4|T(m) + |£ 2 |r(n 2 ). (2.5) 

Since r has been defined as a positively homogeneous function, it is immediate to see 
that (12. 4p for all choices of ni, n 2 , £\ and £ 2 is equivalent to 

r(x + y) < r(x) + r(y) Vx,y. (2.6) 

From (12.61) one easily gets (I2.2p for any tetrahedron. Conversely, if (I2.2p is valid for any 
tetrahedron, then we get (12. 6p by deforming a tetrahedron into a prism, sending one of 
the vertex to infinity. The important inequality is (12.61) (or (12.41) ). 

To summarize, for two distinct phases at equilibrium, the interface free energy is a 
continuous convex function, which is positive and sublinear, that is 

(a) t(x) > x ^ , 

(b) r(tx) = t r(x) Vx and all t > , (2.7) 

(c) r(x + y) < r(x) + r(y) Vx,y. 

By a classical result of Minkowski it is the support function of the convex set W T , that is 



r(x) = sup{(x|y) : y G W T } 



(2.8) 
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We say that an interface perpendicular to n is (thermodynamically) stable if 

t(x + y) < t(x) + r(y) Vx, y linearly independent, such that x + y = n . (2.9) 

In that case ( 12. 3ft and ( 12.51) hold. In general the choice of the normal to the interface 
does not matter, so that we also have r(n) = r(— n). This means that r(n) defines a 
norm on E 3 . 



2.2 Interface free energy and equilibrium shape W T 

We suppose that r is given, verifying properties a), b) and c) of ( 12.71) (we do not assume 
that r(n) = r(— n)). Under these assumptions W T is a convex body, i.e. W T C E3 is a 
bounded, closed convex set with non-empty interior. The point is an interior point of 
W T , because the continuity of r implies that r(n) > a > for all n, ||n|| = 1. 

2.2.1 Equilibrium shape W T and its polar dual W* 

A face F of a convex set K C E 3 is a convex subset F C K, such that 

x , y G iT and -(x + y) G F imply x , y G F . 

Two-dimensional faces are called facets of K and one- dimensional faces are called edges 
of K. An extremal point z of K is a face F = {z}. Equivalently, z is an extremal 
point if it cannot be written as z = Ax+ (1 — A)y with x, y G K and A G (0, 1), except by 
taking z = x = y. The set of extremal points is denoted by extK. Minkowski's theorem 
states that a convex body K is equal to the convex hull of its extremal points, that is 

{k k ^ 

x : A^Xj : Aj > , Aj = 1 , Xj G extK Vz , k arbitrary > . 
i=l i=l ) 

Remark. In dimension two we call "facets" the one-dimensional faces of K C E 2 . 

To study the convex set W T we introduce its polar dual W*. The definition of W* 
is based on the dual relationship between non-zero vectors v and closed half-spaces v* 
containing the origin, 

v* := {x: (v|x) < 1} . 
The polar dual or polar set W* of W T is 

W* := p|{x*: x G W T } = {u : (x|u) < 1 Vx G W T } . 

From the definition it is immediate that W* is a convex set. It is a closed set, since the 
scalar product is continuous, and is an interior point of W*. Hence W* is a convex body. 
Since for x G W T we have (x|u) < 1 for all u G W*, we deduce that W** D W T . Since 
W T is a convex body having the origin as an interior point, we also have W T = W**. 
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Indeed, suppose that z G" W T . Since W T is a closed convex set, there exists a plane 
{x: (x|u) = a} which separates z and W T , i.e. 

W T C {x: (x|u) < a} and (z|u) > a. 

Since is an interior point of W T , we can choose u so that a > 0. Therefore, for all 
y £ W T , (y|u/a) < 1, i.e. u/a G W*. Since (z|u/a) > 1 we conclude that z G" W**, so 
that 

w T = w** . 

There is a simple description of W* in terms of r. If r(u) < 1, it follows from (12. II) that 
u G W 7 ^. Conversely, since r is the support function of W T , there exists for any u a point 
z G W T such that (z|u) = r(u). Therefore, when u G W*, we have r(u) = (z|u) < 1. 
Hence 

H/; = {u : r(u) < 1} and r(x) = min{t > : x/t G W^} . (2.10) 

The statements in (12. lOf) mean that r is the gauge function of W*. We can interpret 
the interface free energy either as the support function of W T , or as the gauge function 
of W*. The boundary dW* of the polar set is simply the level-1 surface of r. Since r 
is positively homogeneous, the level-t surface is obtained from the level-1 surface by a 
dilation of factor t. Since (dW*)* = W** and n* = H(n) for any n G dW*, the boundary 
points of W* give a natural labeling of the support planes of W T (see also section 12.2.21) . 




Polar set and equilibrium shape. 2D Ising model J = 1, (5 = 3. 



Remark. Since W T = W** we can exchange the roles of W T and W*. Let p be the 
support function of W*. This function is the gauge function of the equilibrium shape 
W T , so that 

iy T = {x:p(x)<l}. 
By definition of the support function (see (12.81) ) we have 

p(x) = sup{ ( x|y ) : y G W*} = sup =}► ( x|y ) < p(x)r (y) . 

y^o r(y) 

Example. The Holder norms on E3 are defined by 




(E?=i l x il) 1/P ifl<p<oo, 
max{|xj| , % = 1, 2, 3} if p = 00. 
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Suppose that r(x) = ||x|| p . Then 

W T = {x: ||x|| g < 1} and W* = {x : ||x|| p < 1} , 

where q verifies 1/p + 1/q — 1. 

There is another way of expressing the equilibrium shape using r. For any y we 
introduce the affine functional 

x ^ y (x) := (x|y) - r(y) . 

The half-space H(n) = {xG R 3 : n (x) < 0}. The equilibrium shape can be written 



W T = {x: sup</> y (x) < 0}. 
y 

Since r is sublinear it is easy to verify that 

^ / \ //rv / \ \ f if x G W T 

t (x) := sup((x|y) - r(y)) = <^ . 

y I oo otherwise. 

Thus £/ie Legendre transform r* of the interface free energy is the indicator function of 
the equilibrium shape W T . 

2.2.2 Extremal support planes and stability of interfaces 

We say that dH(n) is an extremal support plane of W T iff n (x) cannot be written 
as 

n (x) = Ci0 ni (x) + C 2 n2 (x) Ci > , C 2 > , 

except by taking n x = t x n and n 2 = t2 n , £1 > and t 2 > 0. Let H(ni) 7^ if(n 2 ), Ci > 
and c 2 > be given. Let n = Ciiii + c 2 n 2 and x e W 7 ,-. Then, by sublinearity of r, 

> 0„(x) = ( x|cini + c 2 n 2 ) - r(cin 1 + c 2 n 2 ) > Ci0 ni (x) + c 2 n2 (x) . 

From this we conclude that dH(n) is extremal iff 

r(n) < Cx7~(ni) + c 2 r(n 2 ) Vn 1; n 2 linearly independent, such that ciiii + c 2 n 2 = n , 

that is, iff the interface perpendicular to n is stable (see (12. 9p ). When the support planes 
dH(n) are parametrized by n e dW*, dH(p) is extremal iff n is an extremal point of 
W*. Indeed, if n = Ani + (1 — A)n 2 is a non-extremal boundary point of W*, then 

1 = r(n) < Ar( ni ) + (1 - A)r(n 2 ) < 1 r(n) = Ar( ni ) + (1 - A)r(n 2 ) , 

so that dH(p) is non-extremal. Conversely, if dH(n) is non-extremal, then there exist 
c\ > 0, c 2 > 0, ni and n 2 so that n = ciiii + c 2 n 2 and r(n) = cir(ni) + c 2 r(n 2 ). Putting 
u = n/r(n), u x = ni/r(ni) and u 2 = n 2 /r(n 2 ), we get 

cir(nx) c 2 r(n 2 ) cir(ni) c 2 r(n 2 ) 

u = — -r^r^ H r^ u 2 and — H = 1 > 

r(n Tin) Tin Tin) 
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so that u G dW* is non-extremal. To summarize, the support plane dH(n) of W T , with 
n G dW* , is extremal iff n is an extremal point of W* . This happens iff the interface 
perpendicular to n is stable. 

We can rewrite (12. ip as 

W T = {x: (x|n) < r(n) , Vn G extW*} . (2.11) 

Only the interfaces perpendicular to n G extW* are stable. For those n, r(n) is well- 
defined and can be measured experimentally. For n G dW*, but not extremal, a flat 
interface perpendicular to n is unstable and is subject to a hill-and- valley formation built 
on the interfaces perpendicular to the m G extW* entering into the extremal decompo- 
sition of n (see |H (1951)| ). 

Remark. For a relation between these deformations of the unstable interfaces and the 
phenomenon of thermal faceting, see [H (1951)] and W (1988) pp. 391-393. It should 



be stressed however that experiments with crystals and interfaces are rarely done at 
equilibrium so that the thermal faceting is basically a non-equilibrium phenomenon. 

Mathematically, the support function of W T is defined for any n, and it can be used 
to extend the definition of r. Equivalently, by taking the convex hull of extW* we can 
define r for all n G dW*. Hence the hypothesis (12. 7p done at the beginning of section 
12.21 are not restrictive. 

Example. There is a tangent plane to W T at x iff there is a unique support plane dH(n) 
containing x; the outward normal to W T at x is well-defined and equal to n. A tangent 
plane is always an extremal support plane. Indeed, suppose that n = Arix + (1 — A)n 2 
with < A < 1 and ni, n 2 G extW*; suppose that dH(n) is the unique support plane at 
x G dW T . Then r(n) = (x|n) and 

1 = r (n) = A(x|m) + (1 - A)(x|n 2 ) < Ar(m) + (1 - A)r(n 2 ) = 1. 

Therefore 



A (r(m) - ( x|m )) +(1 - A)(r(n 2 ) - ( x|n 2 )) = 



>0 



and T(rij) = (x|rij) so that dH{xij) is a support plane at x. Hence n = ri! = n 2 , the 
decomposition of n is trivial and n is an extremal point of W*. We have extremal support 
planes which are not tangent planes when W T has an edge or a corner. 



2.2.3 Subdifferentials of r 

The subdifferential of r at x is the set 

dr{x) := {u G E 3 : r(y) > r(x) + ( u|y - x ) Vy} . 

It is a closed convex set. An element of <9r(x) is a subgradient of r at x. One can 
prove that differentiability of r at x is equivalent to the uniqueness of subgradient; the 
unique subgradient is the gradient of r at x, that is, <9r(x) = {gradr(x)} and 

<9t(x) <9r(x) <9r(x)\ 
dxi ' dx 2 ' dx 3 J 



gradr(x) 
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For x G dW*, gradr(x) is normal to the tangent plane of the polar set W* at x. Since r 
is sublinear, 

0r(x) = 0r(tx) Vt>0. 

From the definition of the sub differential of r at x = one gets another characterization 
of W T (r(0) = 0) 

W T = 0r (0) . 

For any n^O, the subdifferential 0r(n) is the subset of the boundary of the equilibrium 
shape W T , which is given by the intersection ofW T and the support plane dH(n), that is 

0r(n) = W T n dH(n) when n^O. 

Moreover, v G 0r(n) iff <9v* is a support plane of W* at n/r(n). 

Indeed, suppose that v G 0r(n). Then 

r(u) > r(n) + ( v|u — n ) Vu. 

Choosing u = 0, we get (v|n) > r(n). Choosing u = 2n and using the homogeneity of 
r, we get r(n) > ( v|n). Therefore r(n) = ( v|n) and v G dH(n). But this also implies 

r(u) > r(n) + ( v|u — n ) = ( v|u ) Vu, 

so that v G W T . Conversely, if r(n) = (v|n) and r(u) > (v|u) for all u, then 

r ( u ) > ( v l u ) = r ( n ) + ( v l u ~~ n ) Vu. 

For the second part notice that if v G 0r(n), then r(n) = (v|n), so that 0v* is a 
support plane of W* at n/r(n). Conversely, if 0v* is a support plane of W* at n/r(n), 
then v G dH(n) and ( v|u ) < 1 for all u G W 7 "^, so that v G W** = W^. 

In summary, a corner of W 7 "^ at x corresponds to a facet of W T perpendicular to x 
since the subdifferential at x is two-dimensional. Edges of W* correspond to edges of W T . 
A facet of W* corresponds to a corner of W T since for each x of the facet of W* gradr(x) 
is the same and each point of the facet labels a different support plane at gradr(x). For 
the other points x G dW* the gradient exists and (gradr(x)|x) = 1; at gradr(x) G W T 
there is a unique (extremal) support plane perpendicular to x G W*. 

Example. For the Ising model, at zero temperature, W T = {x: < 1} and W* = 

{x : ||x||i < 1}. In dimension d there are only 2d extremal points for W* corresponding to 
the 2d facets of W T . All other interfaces are unstable. For d — 2, at non-zero temperature, 
the corners of W* are smoothed out, so that W T has no facet. Moreover, W* is strictly 
convex, so that W T has no corner and all points of dW* are extremal points, hence all 
interfaces are stable. The disappearing of facets is called roughening transition; here 
the temperature of the roughening transition is zero. In the rough phase, i.e. above 
zero temperature, r is differentiate. This situation is generic for two-dimensional cases 
at non-zero temperature. For d — 3, if the temperature is strictly positive and strictly 
smaller than the two-dimensional critical temperature, inequalities fll.7p and fll.Sp imply 
that there are facets for W T corresponding to corners for W*. Notice that by symmetry, 
instead of taking n{6) = (0, — sin^, cos 9) as in subsection 11.2.21 we may also choose 
n(0) = (— sin^, 0, cos6*) or n(9) = (— sin^, cos6*, 0). At a higher temperature the system 
undergoes a roughening transition with the disappearing of facets, but this result has not 
yet been mathematically established for the Ising model. 
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2.2.4 Radius of curvature and stability of interfaces 

In this subsection we assume that r is strictly convex and smooth (except at the origin). 
At each point of dW T there is a well-defined normal n and W T D dH(n) = {gradr(n)}. 
We consider the two-dimensional case W T C E 2 . We set for a G [0, 2tt), 

n(a) := (cos a, sin a) and m(a) := (— sin a, cos a) . 

The vectors n and m are orthonormal and we can decompose gradr(n) as 

gradr(n) = ( gradr(n)|n }n + ( gradr(n)|m }m = r(n)n + ( gradr(n)|m )m , 

since r(n) = ( gradr(n)|n ). We also have 

^-r(n(a)) = (gradr(n(a))|m(a) ) . 
Therefore we can parametrize dW T by 

a i— > gradr(n(a)) = r(n(a)) n(a) + — r(n(a)) m(ai) , a G [0, 2n) . 

da 



From this expression one obtains 

d ( d 2 

— gradr(n(a)) = f r(n(a)) + — r(n(a)) ) m(a) 



p( a ) = r(n(a)) + — r(n(a)) . (2.12) 



and the radius of curvature at gradr(n(a)), which is given by 

da 2 

Strict positivity of p(a) means absence of corner at gradr(n(a)). There is an interesting 
inequality, noticed by Ioffe, which strengthens the stability condition (12.91) . If the radius 
of curvature p(a) is bounded below uniformly, 

inf p(a) = x > , 

a 

then 

r(x) + r(y) - r(x + y) > x (||x|| + ||y || - ||x + y||) Vx, y . (2.13) 

The constant x is the best possible constant. 

In the three-dimensional case, at gradr(n) G dW T , n is normal to dW T . To study 
the curvature at this point one can slice W T by planes containing n. This reduces the 
problem to two-dimensional situations. If the radius of curvature in each slice is bounded 
below uniformly by x > 0, then we get a stronger version of stability inequality (12.31) : for 
any tetrahedron T(A , Ai, A 3 , A 3 ) 

IAxItM + |A 2 |r(n 2 ) + |A 3 |r(n 3 ) - |A |r(n) > xd^il + |A 2 | + |A 3 | - |A |) . 

Indeed, for x = |Ai|ni, y = |A 2 |n 2 and z = |A 3 |n 3 

r(x) + r(y + z) - r(x + y + z) > x(|| x ll + ||y + z ll - ll x + y + z ll) 

and 

r(y) + r(z) - r(y + z) > x (||y|| + ||z|| - ||y + z||) . 
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2.2.5 Isoperimetric inequality 

Gibbs (1878) and Curie (1885) studied a special case of the following variational problem 
concerning the equilibrium shape of a crystal, and Wulff (1901) gave the geometrical 
interpretation of the solution showing that the equilibrium shape is obtained by the 
construction leading to (12. ip . For this reason the equilibrium shape W T is often called 
Wulff crystal. 

Suppose that V C R d , d = 3 or 2, and that n(s) denotes the outward unit normal to 
its boundary dV at s. The surface free energy of this set is given by the surface integral 



(dTi^ 1 is the (d — l)-Hausdorff measure in R d .) The problem considered by Wulff was 
to determine the set V, which minimizes the functional J- T (dV) among a class of subsets 
with fixed volume. We state an isoperimetric inequality which gives the solution to this 
problem. Roughly speaking, the subsets which can be considered are those subsets for 
which the functional T T is well-defined. Denote by \C\ the volume of the subset C C R d . 



Equality holds if and only iff V = W T up to dilation and translation. 
2.3 Summary 

The subject of this article is the definition of the interface free energy r and its thermo- 
dynamical properties once the interatomic interactions of the system are given. Provided 
that one can construct a macroscopic state with a planar interface perpendicular to n, 
one can use formula (11.31) to obtain r(n). The fundamental property of the interface free 
energy is that it is a convex function. The interface free energy can be measured experi- 
mentally at equilibrium only for the interfaces which are thermodynamically stable. By 
convention the physical value of the interface free energy r(n) is given for a unit vector n. 
But, using the extension of r as an homogeneous function, this function can be interpreted 
either as the support function of the equilibrium shape W T = {x: (x|n) < r(n) , Vn}, 
or as the gauge function of W* = {x: r(x) < 1}. Stable interfaces are labeled by the 
extremal points of W*. 

The interface free energy is considered here only from the macroscopic viewpoint; for 
that reason the macroscopic states of are the relevant states. These states should not 
be confused with the infinite-volume Gibbs measures, which describe the states of the 
system at the microscopic scale, when one chooses as reference length the lattice spacing. 
For example, for the two-dimensional Ising model one has macroscopic states with stable 
interfaces for all directions n, while there are only translation invariant infinite-volume 
Gibbs measures. The study of interfaces at the microscopic scale, when the interfaces 
are stochastic geometrical objects, is of course an important, related but different topic 
(see references for section [TJ). When one studies interfaces at the microscopic scale, it is 
often natural and more convenient to replace the definition of the interface free energy by 
another definition, which is based on specific microscopic properties of the system. One 
essential point is to prove that the definition used is the same as the definition (11.31) . 




Then 



T T {dV) >d\W n 



l/d\y\(d-l)/d 



(2.14) 
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3 Bibliographical Notes 

Section [Q |H (1953)| , [RW (1984)| , [W (1988)] , |Z (1988)] are rev iews of physic s on 
interfaces and equilibrium shapes of crystals. The reviews |A (1986)| and |P (1983)| are 
more specifically concentrated on the interface free energy; they are reviews of mathe- 
matical physics. The excellent review [A (198 6) | also contains a lot of information on 
related subjects. Some results concerning the study of interfaces at the microscopic scale 
are discussed in |P (1983)] . The book Pr (2009")] contains a lot of material related to the 



interface free energy, which is not presented here. 

The problem of the spatial distribution of the coexisting phases, starting with the 
monograph [DKS (1992)] , and the papers [ACC (1990)] and |P (1991)1 , has been one of 



the major research themes in mathematical statistical mechanics during the last decade 
of the XXth century, for which deep mathematical results have been obtained for several 
models. The interface free energy (as defined in section 11.2.11) enters in an essential way 
in the formulation of the results, and interfaces are studied at the microscopic scale. See 
the review paper [BIV (2000)] , the mathematical monog raph |C (2004)| and [Pr (2009)] . 



Section 11.21 In relation with the topic of the article, chapter 6 of |G (1999)] is a good 
introduction to the Ising model, but contains few recent results. An up-to-date good 
reference is |V (2009)| . 

Section 11.2.11 The macroscopic limit for the two-dimensional Ising model is discussed 
in |PV (1999)] . Using the methods exposed in BIV (200~0~)] analogous results can be 



obtain for the three-dimensional case. The up-to-date reference concerning proofs of 
existence and convexity of surface tension for ferromagnetic models is MMR (1992)] . 



For more general models see BP (2003)] . The role of the complete wetting of the walls 



in the definition ( 11.4p has been explicitly emphasized in PV (1999)] . The physics of a 



situation where the interface is partially pinned to the wall had been studied for the 
two-dimensional Ising model by Abraham and Ko already i n |AK (1989)| . Mathematical 
results on wetting phenomenon for Ising systems are in [FP (1987)] and |PV (1996)] . 
Wall free energies are treated in |FC (1977)] . 



Comparisons of several definitions of the interface free energy are carefully discussed 
in |A (1986)] and references can be found there. The computation of Onsager is in 
his famous paper where he solved two-dimensional Ising model at zero magnetic field 
|Q (1944)] ; the computation of Abraham and Martin-L6f is in |AM (1973)| . The com- 



mutation of the full interface free energy is detailed in the treatise of McCoy and Wu 
MW (1973)| . A reference where the results of the computation can be easily found is 



RW (1981)1 - 



Section E221 Inequality (EL5J) is pr oven in [BLP (1980)] and QUE) in |LP (1981)] . Step 



free energy is discussed in |A (1986)] and in [BEF (1986)] . Inequality (O) is proven m 



BEF (1986)1 and inequality flL8]) in |BFL (1982)] T 



Section 12.11 The basic reference is |H (1951)] . 

Section 12.21 The basic reference is |H (1951)] , but Herring uses the surface tension 
plot, which is the set of points 



{x G E 3 : x = r(n)n, ||n|| = 1} , 
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for studying r. This is the standard way of presenting r in physics. One gets the 
surface tension plot from dW* by an inversion on the unit sphere (or the unit circle 
in dimension 2). Affine parts of dW* become spherical parts, or circular parts, of the 
surface tension plot. Using the surface tension plot the fundamental convexity property 
of r is hidden (this property is not stated explicitly in |H (1951)] ); one also looses the 
natural labeling of the stable interfaces by the extremal points of W* of subsection I2.2.2L 
The duality between non-zero x and half-spaces x* is the natural language for discussing 
the geometry of the equilibrium shape. The use of the polar set W* as an alternative to 
the surface tension plot is mentioned in MMR (1992)]. The surface tension plot of the 



two-dimensional Ising model is analyzed in [RW (1981)| . 

Inequality (12.131) appears in [I (1994)| in a slightly different form; the formulation 
given is taken from |PV (1999)] . 

The mathematics of subsections 12.2.11 to 12.2.41 is classical and treated thoroughly in 
M (1911)] . An accessible modern reference is chapter 1 of |S (19 93) | . 



The isoperimetric inequality is a classical topic in analysis. There are many papers on 
inequality (j2TTH . e.g. [D (1944)] (one of the early proofs), |T (1978)] and |F (1991)1 - The 
proof of (I2.14p in these papers is based on Brunn-Minkowski inequality. In dimension 
two a different approach is possible DP (1992)] . 
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